Orbital differentiation and the role of orbital ordering in the magnetic state of Fe 

superconductors 
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We analyze the metallic (tt, 0) antiferromagnetic state of a five-orbital model for iron supercon- 
ductors. We find that with increasing interactions the system does not evolve trivially from the pure 
itinerant to the pure localized regime. Instead we find a region with a strong orbital differentiation 
between xy and yz, which are half-filled gapped states at the Fermi level, and itinerant zx, 3z^ —r^ 
and — y^ ■ Orbital ordering between yz and zx orbitals arises as a consequence of the interplay 
of the exchange energy in the antiferromagnetic x direction due to yz localization and the kinetic 
energy gained by the itinerant orbitals along the ferromagnetic y direction with an overall domi- 
nance of the kinetic energy gain. We argue that iron superconductors are close to the boundary 
between the itinerant and the orbital differentiated regimes and that it could be possible to cross 
this boundary with doping. 

PACS numbers: 75.10.Jm, 75.10.Lp, 75.30.Ds 



There is a strong interrelation between the orbital de- 
gree of freedom, the magnetization, and the lattice struc- 
ture in the Fe-superconductors. Unveiling the nature of 
these connections would define the landscape from which 
superconductivity emerges in these materials. One of 
the important issues is the determination of the strength 
of the interactions. Unlike the cuprates, which are an- 
tiferromagnetic Mott insulators when undoped, the Fc- 
superconductors are antiferromagnetic metals, highlight- 
ing the relevance of the itinerancy of the conduction elec- 
trons. Undoped materials have to accomodate six elec- 
trons in the five Fe-d orbitals, with an average filling of 
1.2, close to the one of doped Mott insulators The 
itinerantiri versus localized^i^ origin of the magnetiza- 
tion has been discussed since the discovery of supercon- 
ductivity in these systems. 

On the other hand there is an increasing evidence for 
orbital differentiation and a possible coexistence of itin- 
erant and localized electrons.—"— Angle Resolved Pho- 
toemission Spectroscopy (ARPES) measurements report 
different renormalization values for the various bands 
close to the Fermi energy depending on their orbital 
character j^^^ Similar qualitative conclusions may be in- 
ferred from Dynamical Mean Field Theory (DMFT) and 
slave-spin calculationsjii^"— 

The possible role of orbital ordering in the mag- 
netism is of present interest. The current debate is fo- 
cused on whether it is the leading instability driving the 
magnetismi^ or it appears as a consequence of the mag- 
netic ordering,ii^— as well as its possible relation to the 
observed anisotropic properties J^i^^ri^ In particular, the 
resistivity in the (tt, 0) antiferromagnetic state was mea- 
sured to be larger in the ferromagnetic y-direction than 
in the antiferromagnetic x-direction^^i^ with a change in 
sign upon hole dopingi^ 

In order to shed hght on the role of the different or- 
bitals on the magnetic state of Fe-superconductors, we 
analyze the metallic (tt, 0) antiferromagnetic state as a 
function of the interactions treated within mean-field. 



Close to the non-magnetic phase boundary, electrons are 
itinerant. An S = 2 state compatible with a local- 
ized Ji — J2 description is found deep in the insulating 
regimci^ With increasing interactions the system does 
not evolve trivially from the pure itinerant to the pure lo- 
calized regime. Instead we find a region with a strong or- 
bital differentiation, see Fig.[T] In this region, zx, Sz^—r^ 
and — are itinerant while xy and yz are half-filled 
and have a gap at the Fermi level. These gapped states 
are reminiscent of the localized electrons discussed in the 
literatureiSr— and for simplicity we refer to them as 
localized. At large values of Hund's coupling the itiner- 
ant electrons are also gapped at half-filling while keeping 
a finite density at the Fermi level. We uncover the dif- 
ferent role that orbitals play in the stabilization of the 
orbital ordering and consequently of the (tt, 0) antifer- 
romagnetic state. We find that while superexchange>iS 
between the localized orbitals contributes to generate or- 
bital ordering between yz and zx, it is necessary to invoke 
the kinetic energy gain of the itinerant electrons along the 
ferromagnetic direction to describe the observed features. 
We analyze this result in connection with the resistivity 
anisotropyi^^i^ We argue that iron pnictides are close to 
the boundary between itinerant and strong orbital differ- 
entiation regimes and that it could be possible to cross 
this boundary with doping. 

We consider an interacting two-dimensional five- 
orbital model for the FeAs layer, as described in Rcf. [20| . 
The Fc orbitals arc defined within the one-iron unit cell 
and hence x and y arc given by the Fe-Fc nearest neighbor 
directions. Only local interactions are included. Interac- 
tions with rotational symmetry can be expressed in terms 
of only two parameters: the intra-orbital Hubbard U and 
the Hrmd's coupling We focus on the metallic (tt, 0) 

antiferromagnetic state and study the phase diagram as 
a function of U and Jh/U with interactions treated at 
the Hartree-Fock leveli^ 

For the tight-binding we use the model described 
in Ref. [33 and obtained within the Slatcr-Kostcr 
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FIG. 1: (Color online) This figure summarizes the main re- 
sults of this work. The colors distinguish the diff'erent (tt, 0) 
magnetic phases of the undoped (6 electrons in 5 d orbitals) 
system as a function of the local interaction parameters U and 
Jh^ The grey area is the non-magnetic region. The blue and 
red areas are magnetic with a high moment (parallel orbital 
moments) and a low moment (antiparallel orbital moments) 
state respectively. The white solid lines on the right sepa- 
rate the metallic {U < 3) from the insulating [U > 3) state. 
We have analyzed the orbital differentiation within the blue 
metallic area. We distinguish two different regions that we 
label, with increasing U , as itinerant and strong orbital dif- 
ferentiation. The strong orbital differentiation region can be 
further splitted by the opening of a gap at half-filling, see text 
for discussion. 



formalism^^ that takes into account the symmetry of 
the orbitals and the lattice. In this model the tight- 
binding parameters are analytic functions of the angle a 
formed by the Fe-As bond and the Fe plane. The result- 
ing bands, their orbital compositions, the Fermi surface, 
and the modifications induced by a are all consistent 
with ab-initio calculations!^ This allows to straightfor- 
wardly explore the effect of the Fe-plane lattice structure 
on the electronic properties. Except when specifically 
stated, the results are obtained for the undoped case 
with six electrons per iron n = 6, and regular tetrahe- 
dra a — 35.3°. 

The main results of our analysis are summarized in 
Fig. [1] on top of the (tt, 0) magnetic phase diagram pre- 
viously reported in Ref. [20|- The system becomes insu- 
lating on the right of the solid white lines. In the grey 
region the system is not magnetic and the red area cor- 
responds to a low moment state, in which Hund's rule 
is violated.— 1^"— The blue high moment state fulfills 
Hund's rule. Deep in the insulating regime, this state has 
a spin S = 2 with filled — and half-filled xy, yz, zx, 
and 3z^ — [23]. In the region of larger U , an increasing 
Jh /U leads to metallicity.— We focus here on the metal- 
lic blue state, on the left of the solid line. Decreasing 
the value of U we find different regions which differ on 



FIG. 2: (Color online) (a) Total density of states for Jh /U — 
0.25 and [/ = 2 eV (red) and U = 2.5 eV (black). The Fermi 
energy is at ixi = 0. For U — 2.5 eV, a gap has opened around 
uj ^ —1.5 eV, which corresponds to half-filling as illustrated in 
(b), where the integrated density of states (filling) is shown. 
This opening of a gap at half-filling in the total density of 
states characterizes the region on the right of the dotted curve 
within the strong orbital differentiation area in Fig. [T] In 
contrast, there is no gap for U = 2 eV. (c) Size of the gap at 
half-filling as a function of U for n = 6. (d) Gap as a function 
of doping n for the two different values of U. The gap has a 
maximum at ?i = 5. 



their electronic structure and the related orbital differ- 
entiation. The regions have been labelled strong orbital 
differentiation and itinerant. The nature of the different 
regions can be inferred from the analysis of the density 
of states, magnetization and orbital filling curves. 

We first focus on the .strong orbital differentiation re- 
gion of the phase diagram. Fig. [2] represents the total 
density of states in two points of this region on both sides 
of the dotted curve in Fig. [TJ [/ = 2 eV and U = 2.5 eV, 
both with Jh/U — 0.25 and n = 6. In both cases, the 
system is metallic with no gap at the Fermi level. How- 
ever, the two curves are qualitatively different with a 
gap clearly showing below the Fermi energy only for the 
largest value of U. The opening of this gap in the phase 
diagram is signaled with a dotted curve in Fig. [1] In 
Fig. mb) it is shown that the gap opens at an energy 
that corresponds to half-filling (five electrons in the five 
d orbitals). This gap was found in LDA-f-U calculations 
and was associated with a Mott gap.— This association 
is compatible with the trends shown by the gap in our 
Hartree-Fock calculation: It increases upon hole doping 
(decreasing n) reaching a maximum at n = 5 (Fig. [2l^d)). 
Moreover, once this gap opens, its size depends linearly 
and much stronger on U than the Slater splittings at 
other energies (Fig. [^IJc)). 

More can be learned about this .strong orbital differen- 
tiation region by looking at the projection of the density 
of states on the five Fe d orbitals. The integrated density 
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FIG. 3: (Color online) Orbital fillings as a function of energy 
in the strong orbital differentiation regime for (a) U — 2.5 eV 
and (h) U = 2 eV. For U = 2.5 eV all the orbitals show the 
gap at negative energies, yz and xy are gapped at the Fermi 
energy while the other three orbitals are gapped only below 
the Fermi energy. In contrast, for [7 = 2 eV only two orbitals 
{yz and xy) show the gap at half-filling while the other three 
orbitals are itinerant. 



of states (filling) of the orbitals as a function of the energy 
is shown in Fig. [3] for the same two points of the phase 
diagram as in Fig. [5] The first thing we notice is that for 
both values of the interaction two orbitals {xy and yz) 
open their gap at half-filling: their gap is already quite 
large for ?7 = 2 eV and extends up to the Fermi level. 
The presence of such a large gap pinned at half-filling 
suggests a localization of xy and yz. This localization is 
not apparent when looking at the total density of states, 
Fig. [2] The opening of the gap at negative energies in the 
total density of states is then related to the appearance 
of the gap in the other three orbitals zx, — y^ ■, and 
3z^ — r^. However, these three orbitals, unlike xy and 
yz, remain itinerant in this region, with a finite density 
of states at the Fermi energy. 

The orbital differentiation just discussed is concomi- 
tant with orbital ordering: zx tends to be more filled 
and away from half-filling than yz, which is stuck to 
half-filling, yz {zx) has a larger intraorbital exchange 
along the antiferromagnetic a;-direction (ferromagnetic 
y-direction) favoring localization of yz (delocalization 
of zx) in a magnetic state with (tt, 0) ordering.— This 
anisotropic exchange comes from the counter-intuitive 
hoppings relation It^^^y^l > \t\,z^yz\ arising from the com- 
bination of direct Fe-Fe and indirect Fe-Pn-Fe hopping 
amplitudes.— On the other hand, exchange between xy 
and yz {zx) is finite only in the ferromagnetic y-direction 
(antiferromagnetic x-direction) and opposes the localiza- 
tion of yz. As a consequence, as shown in Fig.|4l[b), there 
is more energy gain due to the localization of zx than to 
the localization of yz for regular tetrahedra (a = 35.3°). 
For elongated tetrahedra (a > 35.3°), the localization 



of zx becomes much more advantageous in terms of ex- 
change energy, but this does not affect the sign of the or- 
bital ordering in Fig. |4|a). Moreover, the trend changes 
for smaller values of a, where the localization of yz brings 
an energy gain, but this is not refiected in the magnitude 
of the orbital ordering in Fig.|4l[a). 

A deeper analysis of the orbital dependent hoppings 
as a function of a helps to clarify the situation. The 
hoppings involving yz and zx are anisotropic in the 
plane,— 1^ see Fig. |4l^c) . Of those, for a regular tetrahe- 
dron (a = 35.3°), the largest hoppings in absolute value 
in the x-direction are 2 2, , 2 2, and 

yz,yz^ yz,x'^—y'^^ yz,3z'^—r'^^ 

^%x,xy By symmetry, simply exchanging yz o zx, the 
largest hoppings in the y-direction are f^zx^zx^ ^^.x x^-y^^ 
^zx 3z2_,-2; and tyz,xy From these relations, we see that 
three orbitals {zx, x^ — y^ , and 3^^ — r^) prefer to be 
itinerant to gain kinetic energy in the ferromagnetic y- 
direction. This gain in kinetic energy is maintained when 
the tetrahedra are elongated while the localization of yz 
would become even more unfavourable in terms of ex- 
change energy. In squeezed tetrahedra the hopping be- 
tween zx and 3z^ — along the y-direction strongly de- 
creases and induces a reduction of the orbital ordering. 
Therefore, the orbital ordering, which does not change 
sign within the experimentally relevant values of a, see 
Fig. HJa), arises due to the interplay of the exchange en- 
ergy in the a;-direction due to yz localization and the 
kinetic energy in the y-direction due to the itinerancy of 
zx, x^ — y^, and 'iz^ — r^. Depending on a, these two 
factors cooperate (squeezed tetrahedra) or compete (reg- 
ular or elongated tetrahedra) , with an overall dominance 
of the kinetic energy gain. Therefore, this kinetic energy 
gain is important to stabilize the (tt, 0) magnetic ordering 
in the orbital differentiation regime. 

The orbital differentiation sustained by the gain in ki- 
netic energy in the y-direction survives upon doping, as 
shown in Fig. [5l With decreasing n (hole-doping) zx 
filling decreases fast, changing the sign of the orbital or- 
dering at some point but, remarkably, it is kept away 
from half-filling, together with 3z^ — r^, even at n = 54^ 

So far, we have discussed the nature of the strong or- 
bital dijferenciation region of the phase diagram in Fig.[T] 
If we go on decreasing U , xy and yz lose the localization 
features, getting a finite density of states at the Fermi 
level, and all carriers are itinerant. 

In Fig. [5] (a) the magnetization is depicted as a func- 
tion of U (black curve) while the evolution of the orbital 
fillings and their derivatives with U are shown in (b) and 
(c) respectively. For the values of U with strong orbital 
differentiation the magnetization has the typical concave 
shape, but for smaller values of U, the shape of the curve 
is more complex, xy and yz orbitals are shown to go to 
half-filling at around the same value of the interactions 
at which the magnetization becomes concave. Therefore, 
we estimate the value of the interaction U* at which xy 
and yz localize from the change of curvature in the mag- 
netization, which is concave for U > U*. By fitting this 
concave part, the value of U* is given by the intercept 
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FIG. 4: (Color online) Orbital ordering rizx — ny^, exchange 
energy, and the relevant hopping amplitudes as a function of 
a, the angle formed by the Fe-As bond and the Fe plane. 
a = 35.3° (highlighted by a vertical dotted line) corresponds 
to regular, a > 35.3° to elongated, and a < 35.3° to squeezed 
tetrahedra. (a) Orbital ordering for different values of U. It 
gets supressed with increasing a but it never changes sign 
within the experimentally relevant values of a. (b) Exchange 
energy in the (tt, 0) state due to the localized orbitals and 
calculated assuming localization of xy and yz (black curve) 
and localization of xy and zx (red curve). As the tetrahedra 
elongates, this contribution to the energy would favor a local- 
ization of zx versus localization of yz, but this does not hap- 
pen, see text for discussion, (c) Relevant hopping amplitudes. 
Note that, by symmetry f^^^^ = tyz,yz, *L,;c2-h2 = ^lz,x'i-y^ . 
*L,3z2_^2 = *az,3z2_^2, etc.— 



of the fitting (red-dashed) curve with the a;-axes. U* is 
represented by a dashed line in Fig. [TJ 

The derivative of the magnetization displayed in 
Fig. inja) shows two peaks. The seeond one coincides 
with the appearance of the concave behavior and the ten- 
dency of xy and yz orbitals to half-filling. The first peak 
is associated with a reorganization of the Fermi surface 
in the ordered state: below, a spin-density wave between 
the hole and electron pockets at (0,0) and (7r,0) respec- 
tively characterizes the magnetic state, while above, an- 
other spin-density-wave instability between the electron 
pocket at (0, tt) and the hole pocket at (tt, tt) gaps this 
region of the Brillouin zone. This hole pocket is at (0, 0) 
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FIG. 5: (Color online) Orbital filling as a function of doping 
n for U = 2 eV and Jh/U — 0.25. yz and xy are closer to 
half filling for all dopings. Doping with holes towards n — 5 
takes all orbitals closer to half-filling but the kinetic energy 
gain keeps 3z^ — and zx itinerant. 




2 2.5 1.5 2 

U(eV) U(eV) 



FIG. 6; (Color online) (a) Magnetization (black curve) as 
a function of the Hubbard parameter U for Jh/U = 0.25. 
Its derivative (blue, in arbitrary units) has two peaks. The 
red-dashed curve is a fitting to the concave part of the magne- 
tization which is concomitant with the orbital differentiation. 
The inset shows the values of the U* at which the magnetiza- 
tion changes from a convex to a concave shape, estimated from 
the point at which the fitted curved intercepts the x-axis. For 
values ofU <U* , the exact shape of the curve depends on the 
Fe-As-Fe angle and the tight-binding model, see Fig. [T] (b) 
Orbital fillings versus U and (c) their derivatives. At around 
U* , xy, followed by yz, tends to half-filling. 



in the 2-Fe Brillouin zone. 

Therefore, the detailed behavior of the magnetization 
in the itinerant region with U < U* depends on the de- 
tails of the Fermi surface and the tight-binding model and 
in particular on the presence or absence of a hole Fermi 
pocket at (tTjTt). To illustrate this, we show in Fig. [7]the 
dependence of the magnetization on U for our model^ 
with different values of a, all with a hole Fermi pocket 
at (7r,7r), and for a different tight-binding model^ that 
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FIG. 7: (Color online) Magnetization versus U for Jh/U = 
0.25 for different values of the angle a formed by the Fe- 
As bond with the Fe plane, and for a different tight-binding 
model reported in Ref. [Ij. For sufficiently high values of U, 
the magnetization recovers a typical concave shape while for 
the smaller values of U the curves have different shapes. 

presents no such Fermi pocket. The Fermi surfaces for 
these different cases can be found in Refs. [s^l and jil ). 
respectively. While both the value of U at which mag- 
netism appears and U* depend slightly in the model un- 
der consideration, the description presented here in terms 
of itinerant and strong orbital differentiation regimes is 
valid in all these cases. Moreover the value of the mag- 
netic moment at the crossover is similar for all these mod- 
els. 

In conclusion, we have found that except for small val- 
ues of the interactions close to the non magnetic bound- 
ary, the metallic (tt, 0) antiferromagnetic state satisfying 
Hund's rule is characterized by a strong orbital differenta- 
tion between half-filled xy and yz orbitals, showing a 
large gap at the Fermi level and itinerant zx, — 
and — y^ orbitals away from half-filling and showing a 
finite density of states at the Fermi level. The large gap 
at the Fermi level shown by the half-filled orbitals points 
towards localization. While within the present approach 
such localization cannot be definitively addressed we be- 
lieve that it will appear in more elaborate calculations. 

The larger tendency to localization of xy agrees with 
it being close to half-filling in the non-interacting bands 
and thus very sensitive to interactions, while 3z^ — 
and — have the largest filling and thus more itin- 
erant behavior. This result is consistent with Dynamical 
Mean Field and slave-spin calculations in the paramag- 



netic statciiii^"— Moreover, our calculations uncover a 
strong orbital differentiation between yz and zx, in the 
antiferromagnetic state where the tetragonal symmetry 
is broken. 

In a previous work, wc found orbital ordering to pro- 
mote a larger conductivity in the ferromagnetic direction 
except in a striped region of the phase diagram close to 
the magnetic transition (see Fig. 1(b) in Ref. [23| ) . This 
is consistent with the stabilization of the (tt, 0) magnetic 
state driven by itinerancy in the y-direction (and the con- 
comitant orbital order) in the region of the phase diagram 
with both itinerant and gapped carriers. 

Current estimates for the interactions^ situate the 
iron superconductors close to the boundary between the 
itinerant and the strong orbital differentiation regimes. 
This is in accord with the values of the magnetic moment 
that we obtain in this region, expected in our approxima- 
tion to be similar to those found in LDA calculations.— 
Moreover, the anisotropy of the conductivity found ex- 
perimentally in the 122 compounds j^i^^ which is larger 
in the antiferromagnetic direction, is consistent with the 
system not being very deep in the orbital differentiated 
regime, where the opposite sign of the anisotropy has 
been calculatedi^ The exact position in the phase dia- 
gram could be different among families. 

Finally, the boundary between the itinerant and the 
strong orbital differentiation regimes shifts to lower 
(larger) values of the interaction with hole (electron) 
doping. Therefore, by changing the doping it could be 
possible to cross this boundary and enter a different 
regime. Electron doping promotes itinerancy while hole- 
doping can induce selective orbital localization. This is 
consistent with recent DMFT calculations^ in the non- 
magnetic state which show a selective Mott transition of 
the xy orbital with hole doping. The selective localization 
of this orbital could had been observed already in ARPES 
measurements 1^ A change in the sign of the resistivity 
anisotropy in the magnetic state with hole doping has 
been recently reportedi^ Whether this is caused by hole- 
doping induced orbital differentiation is at present not 
known. 
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